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NOTATION 


AB distance from end of runway to touchdown, m 

AC distance from end of runway to start of takeoff roll, m 

Lrw runway length, m 

£ n direction cosine of nth segment in horizontal plane 

R slant range from aircraft to observer, m 

s n length of nth segment in horizontal plane, m 

T 

U position vector (u v w) , m 

T 

X position vector (x y) , m 

T 

X' position vector (x' y' ) , m 

T 

X" position vector (x" y") , m 

T 

0C position vector (x' y' z') , m 

x^£ final coordinate of nth segment, m 

. x' . initial coordinate of nth segment, m 
ni 

y flight-path angle between velocity vector and horizontal plane, deg 

A translation transformation from x",y !l to x,y coordinate system 

0 change in heading angle between segments (positive clockwise or to the 
right), deg 

A rotation transformation from x",y" to x' ,y' coordinate system 

ip heading angle, angle between the x axis in horizontal plane and projec- 

tion of V in horizontal plane (positive counterclockwise) , deg 

0. transformation from u,v,w to x',y',z' coordinate system 


v 




A CLOSED-FORM SOLUTION FOR NOISE CONTOURS 


Elwood C. Stewart and Thomas M. Carson 
Ames Research Center 


SUMMARY 


In this report is described an analytical approach for generating noise 
contours that overcomes the difficulties of existing programs. Although many 
noise-contour programs are available, they require large computers and a sus- 
tained major effort for use, are complex and slow, and are generally not 
desirable for interactive use for on-line piloted-simulator use. The method 
developed here is valid for arbitrarily complex paths and reveals the impor- 
tance of various factors that influence contour shape and size. The calcula- 
tions are simple enough to be implemented on a small, hand-held programmable 

calculator, and a program for the HP-67 calculator is illustrated. The method 
is fast, simple, and gives the area, the contour, and its extremities for 

arbitrary flight paths for both takeoffs and landings. 


INTRODUCTION 


Noise impact on the community is one of the most serious problems of air- 
craft operations. The importance of this problem has been highlighted by the 
increased public concern with environmental issues. To deal with the problem, 
NASA is engaged in research on technical approaches for noise reduction appli- 
cable to new aircraft. 

A key tool in the study of noise impact has been noise-contour calcula- 
tion programs, a number of which have been developed over the last 10 years 
(e.g., refs. 1-3). These programs and their many modifications have been 
widely used by the Government and by industry. The programs handle very gen- 
eral situations which can include excess ground attenuation, shielding, multi- 
ple flights, and various mixes of aircraft. However, there are many diffi- 
culties in the use of these programs. First of all, the programs are 
generally large and not well documented. Consequently, they require extensive 
study to operate. Second, the programs are difficult to make operational, 
especially if they have been stored for a long time. As a consequence of 
these two difficulties, the successful use of these programs usually requires 
a large and sustained commitment. Third, the programs require a large com- 
puter and it is often difficult to obtain sufficient time on a large computer 
for use in an interactive manner. Further, the programs are slow in operation. 
And finally, the programs are generally not desirable for on-line use with 
piloted simulators because of their size and slowness. For these reasons, 
another method that generates noise contours much more rapidly and simply, 
and one that utilizes minimum computational facilities is desirable. 

The present report is concerned with a closed-form analytical approach 
for calculating noise contours that overcomes the difficulties just discussed. 



As is common with most noise programs, an isotropic noise source is assumed, 
although it is believed the present approach is extendable to an anisotropic 
noise source. The method is valid for arbitrarily complex flight paths. For 
simple paths, the calculations can be accomplished by hand; for complex paths, 
the calculations are best done on a small hand-held calculator with program- 
ming capability, such as the HP-67. The analysis reveals the fundamental 
nature of the contours and the importance of the various factors that influ- 
ence its shape and size. The method is fast, simple, and gives the area, the 
actual contour, and its extremities for arbitrarily complex flight paths for 
both takeoffs and landings. 


THEORY 


General Takeoff Solution 

In this section we will consider the noise impact for an aircraft flying 
a takeoff trajectory. A typical three-dimensional trajectory consisting of 
straight-line segments is illustrated by the plan view in sketch (a) below. 


n+1 




Sketch (a) 


Here, way-point 1 is the start of the takeoff roll located a distance AC 
from the origin, way-point 2 is lift-off, and the general ntn segment is 
between the nth and (n + 1) way points, as indicated. Each segment is char- 
acterized by: (1) an angle measured positive counterclockwise from the 

x axis to the segment or, alternatively, by an angle 0 measuring the incre- 
mental change in heading between segments, positive to the right or clockwise 
so as to correspond with existing nomenclature in contour programs (see p. 5 
of ref. 4); (2) an inclination angle y measured positive above the horizon- 
tal; (3) a segment length s; and (4) a distance R from the aircraft corre- 
sponding to the desired EPNdB noise contour and specified thrust level. 

In the following paragraphs interest will focus on determining the noise 
contour, the contour area, and the extremity of the contour, that is, where 
the end of the contour is located. 


Noise contour- The noise contour will be composed of contributions due to 
each of the various segments, and the various pieces of the contour must be 
determined and suitably pieced together so as to generate the complete con- 
tour. For this purpose, consider the path of the aircraft along the nth 
segment in sketch (b) . The aircraft, when represented by an isotropic noise 
source, will generate a constant noise-level surface in the shape of a 
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cylinder as indicated by the following sketch (b) . Now define the orthogonal 
coordinate system u, v, w such that the u axis is aligned with the nth 
segment and the origin is located at the intersection of the u axis and the 
horizontal plane, and the v axis lies in the horizontal plane as shown in 
sketch (b). Another orthogonal coordinate system x' , y' , z' with the same 
origin is shown in which x' is the projection of the u axis on the hori- 
zontal plane, and the y' axis is coincident with the v axis. 



Sketch (b) 


The equation for a general quadric is 

U T ^U = R 2 (1) 

T 

where U = (u v w) . For a cylinder aligned with the u axis. 


a = 


0 0 0 
0 10 
0 0 1 


( 2 ) 


T 

If we define 9C = (x' y' x' ) , the transformation between two vectors QC and 
U is 


0C= ftU 

where ft is an orthogonal matrix given by 

cos y 0 -sin y” 

ft = 0 1 0 

_sin y 0 cos y 

Thus equation (1), when transformed by equation (3), gives 

= r 2 


(3) 


(4) 


(5) 


3 



where 


ado} = 


. 9 

sxn^ y 


0 

1 


[jsin y cos y 0 


-sin y cos Y 

0 


cos 2 y 


( 6 ) 


Equation (5) is a quadric equation for the three-dimensional noise contour for 
a given noise level in terms of the x' , y' , z' coordinate system. 

The intersection of this surface with the horizontal plane will give the 
noise contour. A plane is represented by: 


m 0C = k 
T 

where m = (mg m 2 m 3 ) , and k is a constant. For the horizontal ground- 
level plane, 

T 1 

m =(00 mg) 

k = 0 

Since S)C is orthogonal to m, it is given by 


(7) 


= (x' y' 0) 


T 


( 8 ) 


(9) 


The contour can now be obtained by determining the intersection of equa- 
tions (5) and (9). Partition the following matrices as shown, 


= 


sin 2 y 

o 

-sin y cos y 








A 

B 

0 

1 

0 


C 

D 

-sin y cos y 

0 

? 

cos y 





"x " 



T 


X’ 

y 

- 


— 


0 

_o_ 


L- — 


3T = 


and substitute these two matrices in equation (5) to obtain 

X ,T AX' = R 2 

where 


( 10 ) 


(ID 


( 12 ) 


X 1 


(13) 
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A = 


sin 2 y 0 


( 14 ) 


Equation 

equation 

ellipse. 

contour. 


(12) represents a conic section. Because the diagonal elements in 
(14) are greater than zero, equation (12) can be recognized as an 
Note that only two parameters, y and R, are required to define the 


It should be pointed out that the form of equation (12) applies to each 
segment, although the two constants may differ on each segment. Further, it 
is clear that each segment contributes only a portion of its conic section to 
the total contour. This portion and the way the various portions are to be 
pieced together will be discussed in a few paragraphs. 

The value of R can be determined from the basic noise-source data which 
are typically given as a function of two variables, the slant range R and 
thrust level T, as illustrated in sketch (c). The value of R is determined 



directly from these data for any desired EPNdB contour and for any value of 
thrust T. This latter value could be estimated, determined from flight data, 
or determined from a control generator as described in reference 4. 


An example of special interest is a flight path consisting of only one 
segment at a constant flight-path angle as indicated in sketch (b) . The con- 
tour in this case is given by equation (12) as 


y' 



(x') 2 sin 2 y 


(15) 


For a given R and y, this equation can be plotted directly as function of 
x' . The maximum dimensions of the contour are obviously 


Tiax 


max 


R 


sm y 
= ±R 


(16) 


Now we will consider how to obtain the total contour by properly combin- 
ing the contributions of the various segments. The situation is illustrated 
in sketch (d) which is an elaboration of sketch (a) . Recall that the contour 
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7 , = 0 7 2 

Si s 2 


Sketch (d) 

in the x' , y' coordinate system is given by equation (12). This equation 
must be transformed to the x, y coordinate system, in common with each of 
the other segments. This is done by rotation and translation transformations. 
As for rotation, the transformation is 

X* = A n X" (17) 


where the vector X" is defined by the coordinate system shown in sketch (d) , 
A n is an orthogonal matrix given by 


and \J> n is defined by 


n 


cos sin 

-sin \jJ n cos 4> n 


n 

*n = -E ®i 

i=l 


(18) 


(19) 


in accordance with discussion following sketch (a). Thus, combining equa- 
tions (12) and (17) we obtain 


X " TA n TA n A n X " = V 


As for translation, the transformation is 

X = X" + A n 


( 20 ) 


( 21 ) 
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where 



Combining equations (20) and (21) we obtain 

(X - A n ) T A n T A n A n (X - A n ) = R^ 


( 22 ) 


(23) 


This is the equation of the contour generated by an aircraft flying along the 
extended nth segment in terms of the basic x, y coordinate system. Note 
that this equation applies to each segment but that the parameters for each 
segment may be different. Also note that in addition to the parameters y and 
R, previously discussed with equation (12), there are now two more quantities, 
A and A, involved. 


Now we will determine some of the quantities that are needed in equa- 
tion (23). First, is the matrix A n T A n A n ; it can be shown from equations (14) 
and (18) that 


A n TA n A n 



(24) 


where 


a n = 1 - cos 2 4i n cos 2 y n 
b n = -sin i|J n cos cos 2 y n 
d n = 1 - sin 2 cos 2 y n 


(25) 


Second, is the vector A which is determined as follows: Let n be the 

n 

segment between the nth and (n + 1) way points, let s n be the segment 
length in the horizontal plane, and let & n be its direction cosine matrix: 


Si 


n 


COS ljJ n 

sin i(J n 


(26) 


where ip n is measured positive counterclockwise from the x" axis to the 
segment as discussed previously following figure 1. Now, from sketch (d) it 
can be shown that for the nth segment 
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II 


n- 1 


A n = 


JL s, - 


'-n 


k=l 


k k tan y 


n-1 

X 

k=i 


S k tan Y k + 


AC 

0 


(27) 


The only other variable in equation (23) is R which has been previously 
discussed. 

Equation (23) can be put in a form more suitable for computation. Since 

” - Ax n 
Y - Ay n 


X - A n = 


(28) 


equation (23) can be shown to reduce to 

a n (x " A xn )2 + 2b n (x “ A xn> <7 " A yn> + d n(y ~ A yn > 2 = r2 ^ 

This is the equation for the contour in a useful computational form. For any 
given choice of x, the corresponding value of y is readily obtained as the 
solution of a quadratic equation in y A yn - It should be noted that for a 
given choice of x, the parameters in equation (29) should correspond to the 
segment closest to x. The appropriate segment can be determined by examining 
the geometry of the situation. Or, it can be implemented more precisely as 
follows. Determine the beginning x n p and the end x n f of the nth segment 
by 

n-l 

*ni = Ac + J2 s k cos K (30) 

k=l 

n 

x^f = AC + XI s v cos K (31) 

k=l k 

Then the nth segment to use must be the one such that 

K ni < x < x nf (32) 

Obviously the x values need not exceed x max , the value of which will be 
discussed in a later section. 

Contour are a- The contour area will be composed of contributions due to 
each of the various segments of the flight path. The first segment, which is 
on the ground (sketch (d)) has an area given by 

ttR^ 2 

A 1 = 2R 1 S 1 + g (33 ) 
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The remaining segments are all above the ground, and it will be necessary 
to consider the area contributed by any individual nth segment. An approxi- 
mation is made here that the total contour of an arbitrary path is the same as 
if the trajectory is "straightened out" so that it lies in one vertical plane. 
As illustrated in sketch (e) , this procedure of straightening out the 



Sketch (e) 

trajectory amounts to assuming the two cross-hatched areas are equal. It is 
clear that this approximation becomes worse as the differences between con- 
secutive ip become large. Even then, however, the error is generally only 
a small portion of the total contour area so that the approximation is 
generally very good. 


Consider now the area contributed by the nth segment illustrated in 
sketch (d) . The area between any two points denoted by lower limit LL and 
upper limit UL is given by 


A^CLL.UL) =21 y' dx' 


J *UL 

LL 


/R n 2 - (x' ) 2 sin 2 Y n + -r 


V 


sin y 


sin 1 


n 


x sin y 

r” 


TJL 


n 


(34) 


-III 


where equation (15) has been utilized. This equation applies to each segment 
that is in the air, and the quantities R n and y n are known for each of these 
segments. The limits, however, need some careful consideration. 


For the special case of only two segments, the area due to the second 
segment (i.e., the straight-line segment at a nonzero flight-path angle) is 
obtained by putting LL = 0 and UL = x^ ax as given by equation (16). The 
result is 


ttR,. 


A ? (0, x max ) 


sin y r 


(35) 


This is the total contour area due to the second segment. 

Consider now the determination of the limits LL and UL in equation (34) 
for the area due to the general nth segment. Let L n be the length of the 
contour under the assumption that the aircraft continues to fly the extension 
of the nth segment. From equation (16), L n is 
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Recalling that the origin of the x' , y' coordinate system in sketch (d) was 
determined by the intersection of the extension of the nth segment and the 
horizontal plane, we can deduce that the beginning and end of the nth seg- 
ment are given, respectively, by 


Xni = 


n-1 

E s k tan Yj 
k=l R * 


tan y n 


(37) 


X 


nf 


*Ai + 


n 


(38) 


Then it is clear that 
LL = ^ni 


x Af if 

L n > *Af 

(contour not ended) \ 

L n if 

L n < <f 

/contour ended and n\ | 



\is the last segment / ) 


In the first case, the end of the contour has not been reached by the end of 
the nth segment, while in the second case it has. Then the area contributed 
by the nth segment is 


^ = x'/Rn 2 - (x') 2 sin 2 Y n + 


V 


sm y 


sin 


-1 


x sm y 


lUL \ 


R 


n 


Jx 


ni 


[C 


(40) 


ni 


where F is the quantity in the brackets. This expression is valid for all 
segments including the last one. However, on the last segment, certain round- 
off errors may occur that cause complications in using equation (40). Thus, 
it is usually better to use an alternative expression that avoids this problem 
for the area due to the last segment: 


= 


•^n 2 

2 sin Y n 



(41) 


Another way to avoid problems due to round-off errors is to slightly reduce 
the upper limit in equation (40) on the last segment. 



The total area of the contour due to the total of N segments is now 
;iven by: 

N 


■ - -EH; 

n—2 nl 


( 42 ) 


Con tour extremity - The end of the contour is of special interest, and it 
;an be readily determined from the equations already given. From the last 
section, on the last segment. 


^hnax ^n 


y ' =0 

-’max 


(43) 


Then using the transformations in equations (17) and (21) 

^max “ \ax + A n 


^n X max + A n 


(44) 


Thus 


Xma X _ cos -sin < ax ^ A xn 

Jmaxj [f ln ^ cos VlL°J |A yn _ 


(45) 


General Landing Solution 

Landing problems can be solved by transforming them to equivalent takeoff 
roblems. The landing situation is shown in the vertical plane in sketch (f). 



Sketch (f) 


here B is the touchdown point and s. is the length of the ground-roll 
egment, at the end of which the aircraft has stopped. Now, make the 
ollowing set of transformations illustrated below by the plan views (sketches 
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HORIZONTAL PLANES 



Sketch (g) 


Sketch (h) 


Sketch (i) 


(g)-(i)). Sketch (g) illustrates a landing situation in the given coordinate 
system. In sketch (h) , the direction of the flight has been reversed and the 

origin and direction of the x axis has been changed. The variables are 

denoted by x L and y^. Also, the labeling of the points A', C' , and D' 

corresponds to this new takeoff problem. From this sketch it can be seen that 


A'C’ = L rw - AB - S j (46) 

and A'C' is the variable now used in the equivalent takeoff problem. Also, 

X L = L RW " X (47) 

Finally, in sketch (i), the entire diagram has been rotated clockwise by it 
and the direction of the y R axis reversed. This has maintained the rela- 
tionship between the flight path and the coordinate axes. Then we have 

y L = -y (48) 

Now the problem can be solved in the form shown in sketch (i) as a take- 
off problem. We need to use the value A'C* in this equivalent takeoff prob- 
lem, and to order the segments in sketch (h) or (i) as a takeoff problem 
starting from C'. Also note that 0 is measured positive clockwise so 
that it can be determined directly from sketch (h) if desired. The takeoff 
equations previously developed apply directly using x R and y R as the varia- 
bles instead of x and y. Finally, the contour in terms of x R and y^ is 
converted to the variables x and y in the original coordinate system by 
equations (47) and (48). The extremity of the contour x Lmax is transformed 
by these same equations. 


SUMMARY OF EQUATIONS 


In this section we merely collect the previous equations and procedures 
for determining contour areas, extremities of the contour, and the contour 
itself for both landings and takeoffs. 
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O Takeoffs 
• Area 

ttR^ 2 

Segment 1 = — 2 — + 2R^s 1 

Segment n for n > 2 

Rn 

n sin y n 


n-l 

> s, tan v. 
k-i k 1 


ni 


tan y 


n 


x nf = x ni + s n 


UL = 


if =* contour not ended 

L n if L n < => contour ended, and n 

is the last segment 


^n 


V 


i/lL, 2 - (x') 2 sin 2 Yn + — : sin 

n n sm Y n 


x sin y 


n 


Rn 


UL 


J x ni 


■H 


UL 


x' - 
nx 


where F is defined as the quantity in the brackets. For the last seg- 
ment, one can use 


uRn 
2 sin 


— -H 

l " Y n U 


X . 

nx 


• Extremity of contour 


For last segment, L n < x^ 


Set 


Tiax 


= L. 


n 


Then 


Xm 


ax 


^Snax X Xmax ^ 
= A + 

y 0 A 

_ max_| L J L 


xn 

yn 
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where 



• Contour equations (n > 2) 

Choose x < Xjnax 
Compute x" = x - A xn 

Then compute y" from 

a n (x") 2 + 2b n (x") (y") + d n (y") 2 = R n 2 

where 

a n = 1 - cos 2 cos 2 y n 
b n = -sin ijJ n cos i|J n cos 2 y n 
d n = 1 - sin 2 \p n cos 2 y n 

Then 

y = y" + A yn 

(For n = 1 and x 1;L < x < x^, y = iRj) 

The choice of x and the parameters for the nth segment must satisfy 

x ni < x < x nf 
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where 


n— i 

*ni - AC + £ s k cos \ 
k=l 

n 

x n£ = AC + £ s v cos K 
k=l k K 

O Landings 

• Contour 

Let x T and y replace x and y in the takeoff equations 
JL Li 

Let A’C’ = L„ T7 - AB - s, 

RW A - 

Choose x 

Compute x L = L Ry - x 

Compute y R from takeoff equations 

Compute y = -y R 

• Extremity of contour 


Xm, 


ax 


L. 


RW 
0 . 


" *L, 


.max 


CALCULATOR IMPLEMENTATION 


The equations for the relevant characteristics of the noise contours can 
be implemented on a small hand-held programmable calculator. In this way the 
necessity to access a large central computer system is avoided. In this sec- 
tion we describe the implementation of the contour calculations for one of the 
common programmable calculators, the Hewlett-Packard HP-67. 

Because of the limitation of programming memory in the HP-67, the pro- 
gram was accomplished in two parts: the contour-area and the contour- 

plotting programs. 

The area program calculates and sums the area of the noise contour due to 
each segment and also calculates some quantities necessary for determining the 
noise contour. The contour-plotting program determines the y values of the 
contour associated with an input value x. Instructions for using each of 
these programs are given in figures 1 and 2 for takeoffs and landings. Some 
discussion of the combined use of the two programs is required because there 


15 



are two possible options. The first option, illustrated in figure 3(a), is 
used when the number of segments in the flight path is four or less. The area 
program must be run first because it calculates quantities that are required 
by the contour-plotting program. Next, the contour-plotting program is loaded 
and the contour is plotted. The second option, illustrated in figure 3(b), is 
the more general case. It must be used when the flight path consists of more 
than four segments; it can also be used when the user wants to postpone plot- 
ting the contour to a later time regardless of the number of segments. With 
this option, after the area program calculations have been completed for the 
first four segments, the area and the intermediate quantities are written on a 
magnetic card. The user then goes back to repeat the procedure with more seg- 
ments until the complete flight path is finished. Later, when it is desired 
to plot the contour, the contour-plotting program and the intermediate quanti- 
ties previously computed are loaded in the calculator and the contour is 
plotted. 

For those interested in the details of the program, a flow chart and a 
listing of the area program are given in figures 4 and 5, respectively; fig- 
ures 6 and 7 provide the same information for the plotting program. 


EXAMPLE CONTOURS 


In this section several example applications of the preceding noise- 
contour method are given and compared with the exact results obtained from a 
large contour program (refs. 1-3) using an IBM-360 computer. 

All of the examples to be discussed are based on the same aircraft for 
reasons of simplicity. We will use a typical current conventional jet trans- 
port aircraft that was used in the study of reference 4. The noise character- 
istics for this aircraft are reproduced in figure 8; linear interpolation is 
used for intermediate thrust values. Six examples have been chosen so as to 
exercise many of the possibilities that can arise. The first three examples 
are takeoffs; they and the defining data are illustrated in figure 9. The 
last three examples are landings as illustrated in figure 10. The variables 
s, y, and 0 were chosen for desired geometry while R was determined from 
the noise curve in figure 8 for the given thrust level. 

The contours for the examples are shown in figures 11 through 15. Ground 
attenuation and shielding have been omitted from these examples. Data from 
two sources are shown. The first source of data was obtained from the IBM-360 
program and is illustrated by the contours. The program is about 100,000 bytes 
(8 bit words) in size, and a substantial amount of time and effort was required 
to obtain these results. The second source of data was obtained from the 
closed-form solution using the HP-67 calculator and is illustrated by the cir- 
cles. Comparison of these two sets of data shows there is excellent agreement. 

The contour areas for the six examples are given in the following table 
in which the results from the large noise contour program (IBM-360) and the 
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Example 

Contour area 

, km 2 (mi 2 ) 

Closed-form 

IBM- 

360 


solution 

program | 

I 

154.93 

(59.82) 

154.9 

(59.8) 

II 

154.93 

(59.82) 

155.4 

(60.0) 

III 

100.49 

(38.80) 

99.5 

(38.4) 

IV 

17.59 

(6.79) 

17.4 

(6.7) 

V 

24.22 

(9.35) 

24.1 

(9.3) 

VI 

154.93 

(59.82) 

155.4 

(60.0) 


closed-form solution of this report are compared. The accuracy of the 
closed-form solution is seen to be within 1% of the IBM-360 program results. 


APPLICATION TO JET STOL FLIGHT PATH 


In this section the utility of the method for examining the effect of the 
fundamental variables on the size of the noise contour is illustrated. For 
this purpose, the problem of the effect of flight-path angle y on the size 
of the noise contour is studied. It can be seen that as the flight-path angle 
y increases, the contour tends to become smaller from geometrical considera- 
tions. On the other hand, the increased thrust required to maintain the 
larger flight-path angle results in greater source noise and hence tends to 
enlarge the contour. These two opposing trends are examined. 

The problem to be considered is shown in sketch (j), where the simple 
straight-line path lies in the vertical plane and the angle y is variable. 


SEGMENT 1 

g 

Sketch (j ) 

We take an example lift-fan aircraft, the characteristics of which were sum- 
marized in reference 4. There are two state variables to this problem, y and 
V, both of which are specified. The control variables are thrust, thrust 
angle, and angle of attack, and they are to be determined so as to achieve the 
specified flight path, subject to constraints on each control. This problem 
was considered in reference 4 where it was shown that the controls comprise a 
redundant set and that the redundancy can be removed by an added condition 
that the noise be minimized. It was also shown that the angle-of-attack con- 
straint required for safety due to gusts was the most critical, so that a 
constraint a < 10° was imposed. The solution given in reference 4, then, 
utilized the Kuhn-Tucker theory to determine the optimal set of controls sub- 
ject to the angle-of-attack constraint. This same method was used for the 
present problem. Because of the varying flight— path angle, however, it was 
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necessary to impose an additional pitch-angle constraint, 8 < 20°, for passen- 
ger comfort. Although the solution is obtained for all the controls, thrust 
is most important because it determines the source noise and hence the slant- 
range parameter R for the selected noise level. 

The results of the calculations for a 90 EPNdB contour are shown in 
figure 16. Here the boundary on the control a is shown by a ii m such that 
both constraints a lim < 10° and e n m < 20° are satisfied. The contour 
area is assumed to be due to the second segment (see sketch (j)), because the 
noise from the first segment is largely confined to the airport boundaries. 
Thus, the noise-contour area is 


A„ = : (51) 

2 2 sm 

The value of R (subscripts omitted for simplicity) depends on the selected 
EPNdB contour (90 EPNdB in this example) and the thrust level. The thrust 
level T, as well as the thrust angle n, have been determined as a function 
of the flight-path angle y using the method of reference 4; the results are 
shown in figure 16. Also shown here is the variation of the parameter R as 
determined from the noise data in figure 8. Finally, the contour area A of 
equation (51) is shown as a function of y in figure 16. 

Several interesting observations regarding these results can be made. It 
can be noted that T does not change greatly because the forces are balanced 
by a combination of small changes in T and n over the domain of y. As a 
consequence, the variation of R is relatively small. Thus, the contour area 
varies roughly inversely with sin y, but increases slightly from this rela- 
tionship at the large values of y. For comparison, the inverse variation 
with sin y is shown by the dotted curve in figure 15. 

From these quantitative results, it is seen that the jet STOL aircraft, 
in utilizing its inherent steeper flight-path capability, produces a reduced 
noise-contour area with increasing flight -path angle y because of two con- 
flicting trends. On the one hand, the geometrical effect of a steeper path 
results in a reduced noise-contour area, and it is the dominant effect. On 
the other hand, the increased thrust required for the steeper paths causes a 
rather minor increase in the contour area because of the rotational degree of 
freedom of the thrust. 


The situation would be somewhat different for a conventional jet trans- 
port airplane where there is no rotation of the thrust. In this case, the 
thrust must be increased linearly with y to achieve that flight path, and 
that results in a nearly linear increase in R with y. Since R enters the 
contour area as the second power, it can be surmised that the contour area 
should increase roughly linearly with y, but with possibly some slight curva- 
ture. This trend has been verified by specific calculations for an example 
conventional jet transport airplane using the method of reference 4. 
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CONCLUDING REMARKS 


I The simplicity of the noise analysis presented here enables one to obtain 

the noise contour, its area, and its extremities for arbitrarily complex flight 
paths for both takeoffs and landings. The method is simple and fast, and 
results can be obtained either by manual calculation or by means of a small 
programmable calculator. The use of large programs operating on large com- 
puters is avoided. The analysis reveals the fundamental nature of the con- 
tours and the various factors that influence its size and shape. 

Ground attenuation and shielding effects have been omitted from the anal- 
ysis because their effects are important only on the initial portion of flight 
and are generally confined to the airport boundary. Analysis shows that these 
effects can be included if desired. 

It might be noted that the single-event contour discussed herein is the 
obvious choice for purposes of minimizing noise impact. The impact of multi- 
ple flights can be handled by an obvious extension of the single-event 
results. For example, a common multiple-exposure measure, noise exposure 
factor (NEF), is related directly by a well-known relationship to the single- 
event results. Consequently, desired NEF contours can readily be generated 
by the present analysis. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, April 6, 1979 
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Figure 1.- User instructions for contour-area program. 
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Figure 5.- Concluded. 
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Figure 6.- Flow chart of contour-plotting program. 
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Figure 7.- Code for contour-plotting program. 
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Figure 8.- Noise characteristics for example conventional jet transport. 
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Figure 9.- Takeoff examples. 
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